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Abstract 

This paper introduces a simple and computationally efRcient algo- 
rithm for conversion formulae between moments and cumulants. The 
algorithm provides just one formula for classical, boolean and free cu- 
mulants. This is realized by using a suitable polynomial representation 
of Abel polynomials. The algorithm relies on the classical umbral cal- 
culus, a symbolic language introduced by Rota and Taylor in [TT], that 
is particularly suited to be implemented by using software for symbolic 
computations. Here we give a MAPLE procedure. Comparisons with ex- 
isting procedures, especially for conversions between moments and free 
cumulants, as well as examples of applications to some well-known dis- 
tributions (classical and free) end the paper. 
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1 Introduction 



Among number sequences connected to random variables (r.v.'s), cumulants 
play a central role, characterizing many r.v.'s found in the practice. More- 
over, due to their properties of additivity and invariance under translation, 
cumulants are not necessarily connected with moments of r.v.'s so that they 
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can be analyzed by using an algebraic point of view. For these reasons, cu- 
mulants have been introduced not only in the context of classical probability 
theory, but also in the boolean and the free context [15]. Since cumu- 
lants linearize convolutions of measures - no matter what framework they 
are referred to - this property allows us quick access to test whether a given 
probability measure is a convolution. The linearity of classical convolutions 
corresponds to independent r.v.'s and the linearity of free convolutions 
allows us to recognize free r.v.'s [15] . The boolean convolution was con- 
structed [H] starting exactly from the notion of partial cumulants, this last 
extensively used in the context of stochastic differential equations. 

In this paper, we propose the classical umbral calculus of Rota and Taylor 
[llj as the most natural syntax for conversions between cumulants (classical, 
boolean and free) and moments, without any sophisticated programming. 

The basic devices of the umbral calculus are essentially two. The first one 
is to represent a unital sequence of numbers by a symbol a, called an umbra, 
that is, to represent the sequence l,ai,a2, ■ ■ ■ by means of the sequence 
1, a, , . . . of powers of a via a linear operator E, resembling the expectation 
operator of r.v.'s. This setting is quite natural in free probability [15j. The 
second device consists in representing a sequence 1, 01,02, ... , by means of 
distinct umbrae, as it happens in probability theory too for independent and 
identically distributed (i.i.d.) r.v.'s. Thanks to these devices above all, the 
umbral calculus was already employed successfully as a syntax for symbolic 
computations in statistics [U [7] by using Maple. Multivariate extensions 
have been given in [5j for cumulant estimators as an alternative to tensor 
methods [2]. 

The algorithm here proposed relies on umbral parametrizations of cumu- 
lants (classical, boolean and free) in terms of moments and viceversa, carried 
out in [8j. In the free context, such a parametrization involves umbral Abel 
polynomials. Here we show that it is sufficient to change a parameter in the 
formula in order to get same conversions in classical and boolean theory. 
So we trace back all the parametrizations to one umbral polynomial. We 
give a suitable expansion of such a polynomial in order to get a very simple 
algorithm for the whole matter. 

The paper is organized as follows. Section 2 is provided for readers 
unaware of the classical umbral calculus. We resume terminology, notation 
and some basic definitions. Section 3 recalls the umbral theory of cumulants 
both classical, boolean and free, as well as their corresponding parametriza- 
tions. In Section 4, we show that all the parametrizations can be recovered 
through only one umbral polynomial for which a fruitful expansion is pro- 
vided. The MAPLE algorithm for this expansion is also given. Comparisons, 
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with procedures available in the literature (see for example [T]), confirm the 
competitiveness of the umbral algorithm. The paper ends with some exam- 
ples of how to use umbral calculus in order to compute classical, boolean 
and free cumulants for some classical probability laws. 

2 The classical umbral calculus 

In the following, we recall terminology, notation and some basic definitions 
of the classical umbral calculus, as introduced by Rota and Taylor in pjj 
and subsequently developed by Di Nardo and Senato in [3] and [1]. We 
also recall those results useful in developing the umbral algorithm, given 
in Section 4. We skip any proof: the reader interested in-depth analysis is 
referred to [3j and [4J. 

Classical umbral calculus is a syntax consisting of the following data: a 
set A = {a, (3, . . .}, called the alphabet, whose elements are named umbrae; 
a commutative integral domain R whose quotient field is of characteristic 
zerc0; a linear functional E : R[A] — > R, called evaluation, such that E[l] = 1 
and 

E[a'P^ E[a']E[(3^] ■ • • E[-f''] 

for any set of distinct umbrae in A and for ,k nonnegative integers 

(the so-called uncorrelation property); an element e G A, called augmentation 
[To] , such that -E[e"'] = (5o,n, for any nonnegative integer n, where 5ij = 1 
if i = j, otherwise being zero; an element u € A, called unity umbra, such 
that = 1, for any nonnegative integer n. 

A sequence oq = 1, ai, 02, . . . in i? is umbrally represented by an umbra 
a when 

E[a']=ai, for i = 0,1,2,.... 

The elements are called moments of the umbra a. This name recalls the 
device, familiar to statisticians, when Oj represents the i-th moment of a r.v. 
X. Similarly, the factorial moments of an umbra a are the elements 

0(0) = 1, a(n) = E[{a)n] for nonnegative integers n 

where {a)n = a{a — 1) • • • (q — n + 1) is the lower factorial. The following 
umbrae play a special role in the umbral calculus. 

Singleton umbra. The singleton umbra x is the umbra such that E[x^ = 1 
and Elx^] = for n = 2, 3, .... Its factorial moments are = (— l)"~"^(n— 

^In many applications R is the field of real or complex numbers. 
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1)!. 

Bell umbra. The Bell umbra (5 is the umbra whose factorial moments are 
all equal to 1, i.e. E[{f3)n] = 1 for n = 0, 1, 2, . . . . Its moments are the Bell 
numbers, that is the n-th coefficient in the Taylor series expansion of the 
function exp(e* — 1). 

An umbral polynomial is a polynomial p G -R[^]. The support of p is the 
set of all umbrae occurring in p. lip and q are two umbral polynomials, then 
p and q are uncorrelated if and only if their supports are disjoint; p and q 
are umhmlly equivalent if and only if 

E[p\ = E[q\, in symbols p q. 

It is possible that two distinct umbrae represent the same sequence of mo- 
ments, in such case they are called similar umbrae. More formally two 
umbrae a and 7 are similar when a" is umbrally equivalent to 7", for all 
n = 0, 1, 2, . . . in symbols 

a = 7-^a"~7" n = 0,1,2,.... 

Given a sequence 1, 01,02,... in R there are infinitely many distinct, and 
thus similar umbrae, representing the sequence. 

Thanks to the notion of similar umbrae, the alphabet A can be extended 
by inserting the so-called auxiliary umbrae, resulting from operations among 
similar umbrae. This leads to construct a saturated umbral calculus, in 
which auxiliary umbrae are handled as elements of the alphabet. 

In the following, we focus the attention on some special auxiliary um- 
brae. We assume {q', a", . . . , a'"} is a set of n uncorrelated umbrae similar 
to an umbra a. 

Dot power. The symbol a'" is an auxiliary umbra denoting the prod- 
uct a' a" ■ ■ ■ a'" . Moments of q'"" can be easily recovered from its defini- 
tion. Indeed, if the umbra a represents the sequence 1,01,02,..., then 
E[[a''^)^] = o^ for all nonnegative integers k and n. 

Dot product. The symbol n.a denotes an auxiliary umbra similar to the 
sum a' + a" + ■ ■ ■ + a'" . So n.a is the umbral counterpart of a sum of 
i.i.d. r.v.'s. Moments of n.a can be expressed using the notions of integer 
partitiorll and dot-power. By using an umbral version of the well-known 

^Recall that a partition of an integer i is a sequence A = (Ai, A2, . . . , At), where Aj are 
weakly decreasing positive integers such that Aj = i. The integers Aj are named 

parts of A. The lenght of A is the number of its parts and will be indicated hj vx. A 
different notation is A = (f^ , 2""^ ,•■■), where rj is the number of parts of A equal to j and 
J"! + '"2 + ■ • • = J'A- Note that r-j is said to be the multiplicity of j. We use the classical 
notation A h i to denote "A is a partition of i" . 



4 



multinomial expansion theorem, we have 



{n.aY ~ ^{n),y^dxax, (1) 

where the sum is over all partitions A = (1^^ , 2^^ ^ _ /) of the integer i, {n)^^ = 
for ux > n, 

^^ = ;:^ (l,)n(20^.... a. = [a']--[(a"n----. (2) 

A feature of the classical umbral calculus is the construction of new 
auxiliary umbrae by suitable symbolic substitutions. For example, in n.a 
replace the integer n by an umbra 7. From ([1]), the new auxiliary umbra 
7. a is such that 

(7.a)* Y,{^Udxax, (3) 

and it is the umbral counterpart of a so-called random sum. By using ([3|), 
we have j3.x = X-P = u as well as {a.xY — In the next section, we will 
see that x-<^ h^-s a different meaning. Moreover we have 

(7./?.ar ^5^7'^dAaA. (4) 

Ahi 

The umbra ^.f3.a is called composition umbra of a and 7 and it is such that 
(7./3).a = 7.(/3.a). The compositional inverse a^~^^ of an umbra a is such 
that a<-^>./3.a = x = a./3.Q;<~^>. 



3 Classical, boolean and free cumulants 

In this section, we recall some results given in [8j about connections between 
moments and cumulants (classical, boolean and free). 

a-cumulant umbra. The umbra x-o^j where x is the singleton umbra, is 
called a-cumulant umbra. Replacing 7 by x and by virtue of ([3]), we have 

ix-ay ~ ^x.^dAOA ^ Y.(-ir'-Hiyx - ly.dxax- (5) 

Ahi Ahi 

Since the second equivalence in ([5]) recalls the well-known expression of 
cumulants in terms of moments of a r.v., it is natural to refer to moments 
of x-ck as cumulants of a. The a-cumulant umbra, usually denoted by k^, 
is deeply studied in [Ij. In particular, if is the a-cumulant umbra, then 
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a = 13. K^. Moreover, by recalling equivalence ([3]) and E[{(3)i\ = 1 for all 
nonnegative integers i, we have 

a' ~ {13. ~ Y.^l3),M^c.)x =^ J] (^Ja. (6) 

AH Ahi 

Theorem 3.1 (Parametrizations). Let be the a-cumulant umbra. For 
i = 1,2, . . . we have 

a* ~ + /?.K<,)*"-^ K*^ ~ a(a - l.a)*"-^. (7) 

The symbol —l.a denotes the inverse of an umbra a, that is the umbra such 
that a H — l.a = e. 

a-boolean cumulant umbra. Let M{t) be the ordinary generating func- 
tion (g.f.) of a r.v. X, that is M{t) = 1 + J2i>i where = We 
have 

M{t) = where H{t) = Y,hit\ 

^ ' i>i 

and hi are called boolean cumulants of X. The umbral theory of boolean 
cumulants has been introduced in p]. In particular, the a-boolean cumulant 
umbra r/^ is such that ElrfJ^ = h^. This umbra corresponds to the compo- 
sition of the umbra a, having moments a* ~ iloj, and the compositional 
inverse of u, having moments n!, in symbols f]^ = u'^~^^ .(].a. By this last 
equivalence, we get 

fil^Y.^-ir'~'iyxldxax. 

Ahi 

The previous equivalence allows us to express boolean cumulants in terms of 
moments. In order to get the inverted expressions, we need of the following 
equivalence a = u.p.f]^, that has been proved in the Boolean Inversion 
Theorem (cfr. |8]). Again we have 

a' ~ ^z^a!c?a (??-)a- 

Ahi 

Observe the analogy between the similarity f/„ = u^^^^./^.q;, and the one 
characterizing the a-cumulant umbra k„ = x.a = u"^^^^ .f3.a. 

Theorem 3.2 (Parametrizations). Letrj^ be the a-boolean cumulant umbra. 
For i = 1,2, . . ., we have 

a' ^%{% + 2.u.p.fiJ-^ t ^ aia - 2.ay-\ (8) 
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a-free cumulant umbra. Let us consider a noncommutative r.v. X, i.e. 
an element of an unital noncommutative algebra A. Suppose (j) : A ^ Cisan 
unital linear functional. The i-th moment of X is the complex number = 
<p{X^) while its g.f. is the formal power series M{t) = 1 + X^j>i rriif. The 
noncrossing (or free) cumulants of X are the coefficients rj of the ordinary 
power series R{t) = 1 + Ei>i''i** such that M{t) = R[tM{t)]. The umbral 
theory of free cumulants has been introduced in [8]. As before, the a-free 
cumulant .^^ has been characterized so that = In particular, by 

using the Lagrange inversion formula [3] , the a -free cumulant umbra is such 
that = a^~^^ , where a^, is the derivative umbra of a, i.e. such 

that ~ na""^. It is quite obvious to observe the difference between 
the previous equivalence and the ones characterizing both the a-classical 
and the a-boolean cumulant umbrae. This is why the computation of free 
cumulants is quite difficult compared with the classical and boolean ones. 
Speicher has found a way to expressing free cumulants {r„}„>i in terms of 
moments {'mn}n>i (and viceversa) by using non-crossing partitions of a set 
|12l I13j . However, the resulting algorithm is quite difficult to implement. 
Bryc in [1] uses a different approach. In the next section, we propose an 
unifying algorithm, which relies on the following parametrizations. 

Theorem 3.3 (Parametrizations). Let be the a-free cumulant umbra. 
For 2 = 1,2,... we have 

a' ~ ^^{^^ + i.:^^y-^ ^\ ~ a{a - i-af'^ . (9) 

The umbral polynomials x{x — n.a)""^ are known as umbral Abel poly- 
nomials. 

4 The umbral algorithm 

The algorithm we propose relies on an efficient expansion of the following 
umbral polynomial 7(7 + 5. 7)'"^ for i = 1,2,... and 5, 7 umbrae. 

Proposition 4.1. If 6, j £ A then 

7(7 + M^"'=^E('^)-m-i^m7/.. (10) 

Proof. By using the binomial expansion and equivalence ([T]), we have 

7(7 + M'-'=^E( ^ll J2 (^Udx7x. (11) 

s=l ^ ^ Ahi-s 
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Suppose we consider the partition /x of the integer i, obtained by adding the 
integer s to the partition A. Then we have 7'*7a = and z^;^ = i/^ — 1. If Cs 
denotes the multiphcity of s in A and denotes the multiphcity of s in 
then nis = + 1. Therefore, we have 

/ i-l \ r _ ^ _ ^ 

where the last equahty comes by multiplying numerator and denominator 
for ijis > 0. Recalling that ^sm^ = i, equivalence (fTOl) follows. □ 

In order to evaluate 7(7+15.7)*"^ via (llOp . we need the factorial moments 
of 5 and the moments of 7. Recall that, if we just have information on 
moments 5* , the factorial moments can be recovered by using the well-known 
change of bases: 

i 

(<5),:^^.(i,A;)5^ 
k=l 

where {s{i, k)} are the Stirling numbers of I kind. In particular equivalence 
pop allows us to give any expression of cumulants (classical, boolean, free) 
in terms of moments and viceversa. 

i) For classical cumulants in terms of moments, due to the latter of ([7]), we 
set 6 = —l.u and 7 = 0. Here we find E[{—l.u)i] = (— l)j = (— 

ii) For moments in terms of classical cumulants, due to the first of ([7]), we 
set (5 = /3 and 7 = k^. Here we know E[{l3)i] = 1. 

in) For boolean cumulants in terms of moments, due to the latter of ([8|), we 
set 6 = —2.U and and 7 = a. Here we find E[(—2.u)i] = (— + 1)!. 

iv) For moments in terms of boolean cumulants, due to the first of ([8]), we 
set 5 = 2.u.p and 7 = Here we have E[{2.u.l3)i] = E[{2.u.p.xy] = 
E[{2.uy] = E[{u + u'Y] = + 

v) For free cumulants in terms of moments, due to the latter of ([9]), we set 
6 = —n.u and j = a. Here we have E[(—n.u)i] = {—n)i. 

vi) Finally, for moments in terms of free cumulants, due to the first of ([9]), 
we set 6 = n.u and 7 = Here we have E[{n.u)i] = (n)j. 

The umbral algorithm in MAPLE is the following: 

y:=combinat ['partition'] (i) : 
umbralg := proc(i,fm,y) 

i! * add(fm[j] * 

mul(((g[x[l]]/x[l] !)-x[2])/x[2] ! , 
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x=convert (y [j] , multiset) ) , 
j=l. .nops(y)) ; 

end: 

In the MAPLE procedure, the factorial moments E[{6)j-i] are referred by the 
vector fm[j] and the moments £^[7*^] are referred by the vector g[k] . 

Table 1 refers to computational times (in seconds) reached by using the 
umbral algorithm and the Bryc's procedure [1], both implemented in MAPLE, 
release 7, when we need free cumulants in terms of moments El- 



i 


MAPLE (umbral) 


MAPLE (Bryc) 


15 


0.015 


0.016 


18 


0.031 


0.062 


21 


0.078 


0.141 


24 


0.172 


0.266 


27 


0.375 


0.703 



Table 1: Comparisons of computational times needed to compute free cumu- 
lants in terms of moments. Tasks performed on Intel (R) Pentium (R),CPU 
3.00 GHz, 512 MB RAM. 

Moments of Wigner semicircle distribution. In free probability, the 
Wigner semicircle distribution is analogous to the Gaussian r.v. in the clas- 
sical probability. Indeed, free cumulants of degree higher than 2 of the 
Wigner semicircle r.v. are zero. The first column in Table 2 shows mo- 
ments of the Wigner semicircle r.v. X, computed by the umbral algorithm. 
They are compared with Catalan numbers Cj (second column), since it is 
well-known that E'fX^*] = Cj and E[X'^^~^'^] = 0. By using equivalence (jH) 
it is straightforward to prove that the umbra corresponding to the Wigner 
semicircle distribution is ?./?.(5, where ? is the umbra whose moments are the 
Catalan numbers and 6 is an umbra having only second moment equal to 1, 
the others being zero. In the next section, we will use again this umbra in 
describing the umbral syntax of a Gaussian r.v. 

Moments of Mar chenko- Past ur distribution. In free probability, the 
Marchenko-Pastur distribution is analogous to the Poisson r.v. in the clas- 
sical probability. Indeed, the free cumulants are all equal to a parameter 
A. The last column in Table 2 shows moments of the Marchenko-Pastur 
distribution computed by the umbral algorithm. 

^The output is in the same form of the one given by Bryc's procedure. 
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i 


Wigner 
r.v. 


Catalan 
numbers 


Marchenko-Pastur r.v. 




2 


1 


2 


A^ + A 




3 





5 


A3 + 3A2 + A 




4 


2 


14 


A^ + 6A3 + 6A2 + a 




6 


5 


132 


A^ + 15A5 + 50A^ + 50A3 + 


15A2 + A 


8 


14 


1430 


A^ + 28A^ + 196A^ + 490A^ 


+ 490A^ + 196A3 + 28A2 + A 



Table 2: Moments of some special free distributions. 



5 Computing cumulants of some known laws through 
the umbral algorithm 

An umbra looks like the framework of a r.v. with no reference to any prob- 
ability space, just looking at moments. The way to recognize the umbra 
corresponding to a r.v. is to characterize the sequence of moments 
When the sequence exists, this can be done by comparing the moment gen- 
erating function (m.f.g.) of a r.v. with the so-called generating function 
(g.f.) of an umbra. The g.f. of an umbra has been defined [4] as the follow- 
ing formal power series 

-j-n 

f{a,t) = 1 + y'an-r- 

n>l 

In the classical umbral calculus, the convergence of the formal power series 
f{a,t) is not relevant. This means that we can define the umbra whose 
moments are the same as the moments of a lognormal r.v., even if this 
r.v. does not admit a m.g.f. So the umbral algorithm allows us to compute 
classical!^, boolean or free cumulants of any r.v. having sequence of moments 
{m„}. Once the umbra corresponding to the r.v. has been characterized, 
we choose 5 = —l.u, — 2.n, —n.u in equivalence pO]) . depending on whether 
we need classical, boolean or free cumulants. 

In the following, we take up again some of the examples given in [1], 
showing how they can be recovered through the umbral algorithm by a suit- 
able characterization of the involved umbrae. 

Poisson r.v. A Poisson r.v. of parameter A is umbrally represented by 

^Recall that, due to their properties of additivity and invariance under translation, the 
cumulants are not necessarily connected with the moments of probability distributions. 
We can define cumulants of any r.v. disregarding the question whether its m.g.f. converges 

[al- 
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the umbra A./3, because f{X.(3,t) = exp[A(exp(t) — 1)], which is the m.g.f. 
of a Poisson r.v. The moments are [3] = E[{X.(3y] = 5'(i, A:)A'^, 

where S{i, k) are the Stirhng numbers of second type. So cumulants of a 
Poisson r.v. can be computed via the umbral algorithm, taking as input 
the sequence of moments i?[(A. /?)*]. If the input is the sequence of factorial 
moments {(— we get classical cumulants; if the input is the sequence 
of factorial moments {(— + 1)!}, we get boolean cumulants; if the input 
is the sequence of factorial moments {{—n)i} we get free cumulants (Tables 
3 and 4 in [5). 

Compound Poisson r.v. A compound Poisson r.v. Sn = Xi + • • • + 

Xn, where {Xi} are i.i.d. r.v.'s and is a Poisson r.v. of parame- 
ter A, is umbrally represented by the umbra X.p.a, because f{\.p.a,t) = 
exp{A[/(Q;, t) — 1]}, which is the m.g.f. of a compound Poisson r.v. The for- 
mal power series f{a,t) corresponds to the m.g.f. of X^. Due to equivalence 
(jl]), the moments are = E[{X.p.ay] = "^x^i ^'^dxax- 
Exponential r.v. An exponential r.v. is umbrally represented by the um- 
bra J because f{u/X,t) = (1 — which is the m.g.f. of an exponential 
r.v. with parameter A > 0. So its moments are = = ^. In order 
to obtain the second column of Table 3 in [l], choose in (fTO]) A = 1 and 
{(— n)j} as factorial moments. 

Uniform r.v. An uniform r.v. on the interval [a, b] has m.g.f. 



pi 6 pt a 

M{t) = , =e*" 

^ ' t{b- a) 



1 



s=t(b—a) 



The umbra, with g.f. is the inverse of the Bernoulli umbra [TT], i.e. 

—l.L. So an uniform r.v. on the interval [a, b] is umbrally represented by the 
umbra a + {b — a)(— l.i). Recalling that l.z.)*] = we get 



a, = E{[a + ib-a){-l.t)f} = Y,([) 



V- / ^ ^ jo - QJ 
a ■' 



j + 1 



In order to obtain the last column of Table 3 in [Ij, choose in (llOp a = 
— 1,6 = 1 and {(— as factorial moments. 

Bernoulli r.v. A Bernoulli r.v. of parameter p G (0, 1) has m.g.f. M(t) = 
q +pe* = 1 +p(e* — 1). The umbra with g.f. M(t) is X-P-(^ (see [1] for more 
details), whose moments are = p. 

Binomial r.v. A binomial r.v. of parameters n, a positive integer, and 
p G (0, 1) is a sum of n i.i.d. Bernoulli r.v.'s. Similarly a binomial r.v. is 
umbrally represented by n.{x-P-f^)- Due to equivalence ([I]), we have = 
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E!{[n-{x-P-(^)Y} = J2xH.i^)i^\'^>^P'^^ ■ order to obtain the second column of 
Table 4 in [Ij, choose {{—n)i} as factorial moments in (jlOp. 
Gaussian r.v. A gaussian r.v. with real parameter fj, and o" > has 
m.g.f. M{t) = exp ^fit + c^y^ • This power series is the g.f. of the umbra 

H + I3.{a6), where 6 is an umbra such that E[6'^] = 1 whereas ElS"^] = for 
positive integers i ^ 2. The following proposition gives the expression of the 
n-th moment of the umbra representing the Gaussian r.v. 

Proposition 5.1. For n=l,2,. . . we have 

[n/2j . 2\k ( ^ 

a. = E[{^. + /3.M)r] = E (t ) (12) 

fc=o ^ ^ 

Proof. By using the binomial expansion, we have 
{^i + ^3.{a5)Y ^ p Q/."-^ gdAM)A. (13) 

Suppose A = (ri,2^2, . . .), then {a6)x = {a6'y''^[a^{6")Y^ ■■■ = a^Sx- On 
the other hand, due to the definition of the umbra (5, we have E[5\\ ^ 
iff the partition A is of type (2''^). Thus, if j is odd, there does not exist 
any partition of this type and so [/?.((t5)]"' ~ 0. Instead, if j is even, say 
j = 2k, there exists a unique partition of type (2^'^) which corresponds to 
r2 = k. For this partition, we have dx = ^^f^l^j^, so that [/^.((T^)] ^'^ ~ ^^2!)'fc°fc! • 
Replacing this last equivalence in ([131), have 

(u + B (a6)V ~ V - ^n-2fc i^fc-K 

+ - 2^^ ^2k\){n - 2k)\^ {2\)^kl 

by which the result follows. □ 

We have proved by umbral tools that moments of a Gaussian r.v. can 
be expressed by using a particular sequence of orthogonal polynomials, the 
Hermite polynomials: 



1 — n V / 



{n)2k n~2k 



k=0 ^ ' 

with X = n and v = —a^. All the properties of Hermite polynomials can 
be easily recovered by using the Gaussian umbra. We skip the details. As 
before, the first column of Table 3 in [1] can be recovered from the umbral 
algorithm by using (jl2p with /i = and a = 1 and choosing {(— ?^)^} as 
factorial moments in ()10p . 
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